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The manipulation of single magnetic molecules may enable new strategies for high-density infor- 
mation storage and quantum-state control. However, progress in these areas depends on developing 
techniques for addressing individual molecules and controlling their spin. Here we report success in 
making electrical contact to individual magnetic N^Ceo molecules and measuring spin excitations 
in their electron tunneling spectra. We verify that the molecules remain magnetic by observing a 
transition as a function of magnetic field which changes the spin quantum number and also the 
existence of nonequilibrium tunneling originating from low-energy excited states. From the tunnel- 
ing spectra, we identify the charge and spin states of the molecule. The measured spectra can be 
reproduced theoretically by accounting for the exchange interaction between the nitrogen spin and 
electron(s) on the Ceo cage. 
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Magnetic molecules provide the opportunity to study, 
at a fundamental level, the origins of ferromagnetism and 
quantum aspects of magnetic dynamics [U [2] . Proposals 
suggest that their spin degrees of freedom could also serve 
as useful qubits for quantum manipulation (3) . Most pre- 
vious experiments on magnetic molecules have examined 
samples consisting of many molecules; however, applica- 
tions and more detailed scientific studies will require the 
ability to address individual molecules. Here we report 
the use of tunneling spectroscopy within single-molecule 
transistors (SMTs), a technique which has been used 
previously to measure vibrational and electronic excita- 
tions m [SI [6] , to achieve electrical contact to individual 
molecules of the spin-3/2 endohedral fullerene N^Ceo 
(Fig. la) and to measure its spin excitations. N^Ceo 
is an attractive model system because of its sim- 
ple spin structure and because of the possibility of con- 
trol experiments employing nonmagnetic Ceo molecules 
[2 [3 [la [III [II [13 IB . N@C6o molecules also have the 
advantage of being stable [15 at the high temperatures 
present during the electromigration process by which our 
molecular-scale junctions are formed |^ |^. Previous 
SMT experiments using less robust Mni2-based magnetic 
molecules detected magnetic signatures in a few samples 
[13 [E], but found that the molecular magnetism was 
usually destroyed during device fabrication [19 . We ob- 
serve that the N^Ceo devices exhibit clear magnetic char- 
acter, in that they have a spin-state transition as a func- 
tion of applied magnetic field. The nature of this tran- 
sition allows us to identify the charge and spin states of 
the molecule in the SMT. The spectra of N^Ceo also ex- 
hibit low energy excited states and signatures of nonequi- 
librium spin excitations predicted for this molecule [20] . 



We associate the existence of a spin transition in N^Ceo 
accessible at laboratory magnetic fields with the scale of 
the exchange interaction between the nitrogen spin and 
electron(s) on the Ceo cage. 

The geometry of our SMT devices is shown in Fig. 
lb. We make the devices by adsorbing the molecules 
onto an initially-continuous Pt wire on top of an oxi- 
dized Al gate electrode and then breaking the wire us- 




FIG. 1: Device geometry, and the spin states of N^Ceo- a. 
Schematic diagram of the N^Ceo molecule, b. Scanning electron 
microscopy image of a sample at room temperature following elec- 
tromigration. The scale bar denotes 200 nm. The circuit schematic 
shows our biasing convention, c. Spin states for N^Ceo and its an- 
ions. The blue arrows represent the Sn =3/2 spin of the N atom. 
The horizontal black lines represent the lowest unoccupied molecu- 
lar orbitals of the Ceo- The possible total spin states of the N^Ceo 
ion are indicated. 
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FIG. 2: Color-scale plots of differential conductance (dl/dV) for 
N@C6o and Ceo single-molecule transistors, a, dl/dV for N^Ceo 
device #1 as a function of source-drain voltage and applied mag- 
netic field at constant value of gate voltage Vg = -891 mV on the 
positive side of the degeneracy point. Yellow triangles indicate a 
change in slope of the ground-state peaks, and gray triangles mark 
the non-crossing levels discussed in the text. The labels designate 
different transitions. In region I, the ground-state transition cor- 
responds to decreasing spin, Sz,q+i — Sz,q < 0, while in region 
II Sz,q+i - Sz,q > 0. b. Same' as a for N^Ceo device #2. The 
maximum of the color scale is 0.15 /h. c, dl/dV as a function 
of V and B at constant Vg on the positive side of the degeneracy 
point for a Ceo single-molecule transistor. The maximum of the 
color scale is 0.02 e'^/h. No change in slope is observed for the 
ground-state peaks of Ceo- 



ing electromigration to form a nm-scale gap. (Addi- 
tional information about fabrication procedures is given 
in the Methods section.) We identify the presence of 
one or more molecules in the gap by the observation of 
gate-dependent Coulomb-blockade transport characteris- 
tics O HI] . Our success rate for obtaining such devices 
was 9/19 for N^Ceo devices and 17/59 for Ceo devices, 
while 0/39 control samples prepared using pure toluene 
instead of a fullerene solution showed Coulomb blockade. 
Six of the N@C6o devices were sufficiently stable for de- 
tailed measurements. 

The most striking difference between the N^Ceo de- 
vices and the Ceo controls is seen in the magnetic field de- 
pendence of the tunneling conductance spectra. In Fig. 2 
we show these spectra as a function of source-drain volt- 
age (V) and magnetic field (B) at dilution-refrigerator 
temperatures (electron temperature ^ 100 mK) for val- 
ues of gate voltage (Vg) slightly more positive than the 
"degeneracy point" - the value of gate voltage for which 
the energy difference between the two accessible charge 
states is zero, so that current can flow even near V = 0. 
Figure 3 shows the conductance for each sample as a 
function of V and Vg near the degeneracy point. (By 
focusing on low-voltage transport near one degeneracy 
point, we ensure that tunneling is occurring through a 
single molecule [5 .) In the N@Ceo devices, as a func- 



tion of increasing B (Fig. 2a,b) the lowest- volt age con- 
ductance peaks at positive and negative bias first move 
apart and then change slope to move back closer together. 
These initial peaks can be associated with the same quan- 
tum transition from the ground state of one charge state 
(charge q) of the molecule to the ground state with charge 
q -\- 1 (one fewer electron) . At the value of B where the 
slopes of the peaks change sign in Fig. 2a, we observe 
level crossings, where a conductance peak that at low B 
is associated with an excited state of charge q-\-l becomes 
the ground state of charge g + 1 at high B. We measured 
ground-state conductance peaks that shift apart then to- 
gether as a function of B on the positive-gate side of the 
degeneracy point in four out of five of the N@Ceo devices 
on which we performed this measurement. (The fifth 
showed no change in slope up to 9 T.) The field values at 
which the crossovers occurred varied from approximately 
1 T to 7 T. 

The level crossing in Fig. 2a shows that the change in 
sign of the ground-state-peak slopes is due to an increase 
in total spin component along the magnetic field direction 
for the lowest-energy charge q -\- 1 state. The slopes of 
the peaks as a function of B imply that at low fields 
Sz,q > 5*2, and at high fields Sz^q < Assuming 
that \Sz,q-\-i — Sz,q\ = 1/2 , the slopes of the conductance 
resonances vs. magnetic field shown in Fig. 2 correspond 
to an electronic ^-factor \g\ = 2.0 ± 0.3. 

In contrast. Fig. 2c shows the tunneling spectrum for 
a Ceo device also taken on the positive-gate side of its 
degeneracy point. Of the five Ceo devices we measured 
as a function of magnetic field, none showed a change 
in the sign of the slopes for the lowest- volt age tunneling 
transition. While the spectra for the Ceo devices often 
contain low-energy excitations in the same range of en- 
ergy as those of N@Ceo, 0.1 meV to several meV, the lack 
of spin transitions suggests that the low-lying excitations 
in Ceo are vibrational satellites [41 [21], rather than 
being electronic or magnetic states. 

Recently, a tunneling transition with a change in slope 
as a function of magnetic field has been observed in a 
Ceo device, for gate voltages in a limited range near the 
middle of one Coulomb-blockade region [M]. Our mea- 
surements showing no such level crossing in Ceo devices 
are not in conflict with ref. [14] because our data for Ceo 
and N@Ceo likely correspond to a different charge state, 
and our measurements were made only for gate voltages 
near the degeneracy point, well away from the middle of 
any Coulomb-blockade region. 

We assert that the spin-state transition in the N@Ceo 
devices is likely associated with two states with differ- 
ent values of total spin, as opposed to two states within 
the same multiplet split by magnetic anisotropy. Because 
any magnetic anisotropy axes could be oriented arbitrar- 
ily relative to the applied magnetic field, if the transi- 
tion were due to states split by magnetic anisotropy one 
should expect, in general, that the slopes of the energy 
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FIG. 3: Color-scale plots of differential conductance (dl/dV) as a 
function of bias voltage and gate voltage at zero applied field (B = 
T). a,b, dl/dV from the same N^Ceo devices #1 and #2 shown 
in Fig. 2a and b, respectively. The labels in a designate different 
transitions, consistent with Fig. 2a. Their energies, determined by 
extrapolating the transition lines to the Coulomb-blockade thresh- 
old, are (± 0.02 meV) Eb-Ea= 0.29 meV, Ec-Ea= 0.48 meV, 
Ed-Ea= 0.93 meV, Ee - Ea = 1.35 meV, E^-Ea = 0.20 meV 
and Ec — Ea = 0.41 meV. The triangles indicate the termination of 
ETE peaks associated with non-equilibrium excitations, c, dl/dV 
for the same Ceo device shown in Fig. 2c. No ETE peaks are 
visible for this nonmagnetic molecule. 

levels as a function of B should show deviations from 
1^1 =2 and exhibit strong deviations from linearity at 
low field [191 llll US] • In our data, the field dependence 
of the transitions is very close to linear (away from local- 
ized kinks associated with level crossings). The lack of 
observable magnetic anisotropy effects is consistent with 
the symmetries of N^Ceo and the weak spin-orbit cou- 
pling of its constituent atoms [24 . 

With this understanding, the existence of a spin-state 
transition at an accessible magnetic field for the charge 
q -\- 1 energy levels of N^Ceo allows us to identify the 
charge and spin states involved in the tunneling. Con- 
sider the possible states shown in Fig. Ic, which depicts 
the nitrogen spin (blue arrow) coupled to electron spins 
(black arrows) occupying the lowest unoccupied molecu- 
lar orbital (LUMO) of the Ceo- We assume that a 3- fold 
degeneracy of the LUMOs of isolated Ceo is broken by 
molecular deformation due to interaction with the elec- 
trodes or the Jahn- Teller effect [25 , and we consider only 
the experimentally accessible neutral and negative charge 
states. Neutral N@Ceo is known to have spin 3/2 be- 
cause the encapsulated nitrogen retains its atomic elec- 
tronic configuration [26 . When electrons are added to 
N@Ceo, they are believed to occupy the LUMO of the 
Ceo molecule [27]. With the addition of one extra elec- 
tron, N@Cgo may therefore have either 6* = 1 or 2, and 
these two multiplets are split in energy by the exchange 
interaction between the N spin and the electron spin. If 
the exchange is antiferromagnetic, the ground state will 



have S' = 1 for low B smd S = 2 for large B. For N@C^o ^ 
the total spin should again be 6* = 3/2 since the splitting 
between the LUMO states is expected to be large com- 
pared to the Hund's rule coupling. Because we observe 
the spin-state transition in the more positive of the two 
charge states involved in tunneling, the 0/1- charge cou- 
ple can be ruled out as the charge states between which 
tunneling occurs. This leads to the identification, for an 
antiferromagnetic exchange interaction, that 1-/2- is the 
couple with the lowest charge magnitude that is capable 
of transitioning from Sq — Sq-^i = +1/2 at low fields to 
Sq — Sq-^i = —1/2 at high fields. For a ferromagnetic ex- 
change interaction, the lowest-magnitude charge couple 
possibly consistent with the data would be 2-/3-, which is 
unlikely because the 3- charge state should have a much 
higher energy than 1- or 2- states. A 2-/3- scenario would 
also require fine tuning such that the exchange splitting 
is approximately equal to the deformation-induced level 
splitting within the N^Cqq LUMO. 

The identification of our measured spectra with the 1- 
/2- charge couple implies that our N@Ceo devices are in 
the 2- charge state at = 0, which may seem to be a 
large magnitude of charge in an unbiased device. How- 
ever, for Ceo in the gas phase, both the 1- and 2- states 
are more stable than the neutral form, with the 2- state 
quite close in energy to the 1- state [28l [29l |30] . The equi- 
librium charge of Ceo is generally found to be either 1- or 
2- when it is adsorbed on noble metal surfaces [SU [32] . 
These previous studies suggest that the 3- charge state 
is likely to have a much higher energy than 1- or 2- in 
an unbiased device, which argues against the possibility 
of a ferromagnetic exchange interaction in our samples 
as this would require a 2-/3- charge couple with a 3- 
charge state at zero gate voltage. In regard to the sign 
of the exchange interaction, an ab-initio calculation has 
predicted a weakly ferromagnetic interaction in N@Cgo ? 
rather than the antiferromagnetic interaction we suggest 
[33 . However, it should be noted that the calculation of 
the exchange interaction is difficult since it is computed 
as a small difference between large total energies. 

The device-to-device variations that we measure in the 
N@Ceo tunneling spectra indicate that the energies of the 
electronic states are perturbed differently by the environ- 
ment in each device. In order to explain the variations 
in field value giving rise to the spin-state transition, the 
strength of the exchange interaction must differ between 
devices, presumably due to different amounts of molec- 
ular deformation and/or local electric fields. If we pa- 
rameterize the exchange interaction in the form Jsg • Sat, 
where Sg is the total electron spin and Sat is the spin on 
the nitrogen atom, then measured transition fields (1 to 
7 T) imply exchange strengths, \J\ , between 0.06 meV 
and 0.4 meV. 

To understand how the measured tunneling spectra 
might relate to the molecular magnetism, we have cal- 
culated the conductance of an N@Ceo transistor (Fig. 
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FIG. 4: Numerical calculation of the tunneling spectrum for N@C eo, with parameters chosen to mimic the data for Device #1. a, 
Calculated dl /dV as a function of V and B at constant Vg on the positive side of the q = -2 / q = -1 degeneracy point. Compare to 
Fig. 2a. The yellow triangles mark the spin-state transition, and the gray triangles mark non-crossing levels, b, Calculated dl /dV as a 
function of V and Vg. The triangle marks the termination of a nonequilibrium ETE peak. Compare to Fig. 3a. c, Energy level diagram 
for a large applied magnetic field showing allowed tunneling transitions which contribute to the low-energy spectrum, for tunneling from 
the charge q = -2 states (blue) to the q = -1 states (red) of N^Ceo- The dotted lines show nonequilibrium transitions. (For B = 0, 
the different ms levels within a given multiplet will become degenerate in the absence of magnetic anisotropy.) d. Allowed tunneling 
transitions contributing to the low-energy spectrum, for tunneling from q = -1 states to q = -2 states. The letters labeling the transitions 
in panels c and d can explain the labeled states in Figures 2a, 3a. The relative energies of the q = -1 and q = -2 states will depend on Vg. 



4a, b). As a minimal model to account for the low-energy 
transitions, we consider the energy level diagrams shown 
in Fig. 4c, d, which include the three multiplets, {q = -2^ 
S = 3/2), {q = -1, 5 = 1), and {q = -1, 5' = 2), we have 
introduced previously and also the lowest-energy excited- 
state multiplet associated with each of these states. The 
excited multiplets could correspond to the added excita- 
tion of a long-lived vibrational mode or energy shifts due 
to charging in nearby molecules. We assume that the 
molecule has negligible magnetic anisotropy, so that all 
of the magnetic effects we consider are due to exchange- 
induced level splittings. This is in contrast to previ- 
ous calculations for Mni2-type molecules where the low- 



energy spectra were assumed to be dominated by large 
magnetic anisotropy [18, 19 . We use the rate-equation 
approach [20 to calculate the differential conductance. 
Parameter values, listed in the Methods section, are se- 
lected to approximate the data in Fig. 2a and 3a. 

There are a sufficiently large number of free parame- 
ters that we do not seek or claim quantitative agreement. 
Depending on assumptions about the value of J and vi- 
brational parameters, different numbers of excited-state 
tunneling transitions are possible. Nevertheless the cal- 
culation can reproduce many of the qualitative features 
seen in Fig. 2 and 3. First, the N^Ceo spectra shown 
in Fig. 3a, b and the theoretical results shown in Fig. 4b 
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both exhibit a pattern of nonequihbrium excitations that 
was predicted [20^ for spin excitations in N^Ceo- If one 
extrapolates the conductance peaks marked by yehow tri- 
angles in Fig. 3a and b toward zero voltage, they do 
not extend all the way to the line corresponding to the 
ground-state tunneling threshold, but instead terminate 
when they intersect a conductance peak corresponding to 
an excited-state transition. We observed this type of fea- 
ture in four of six N^Ceo devices, but in none of thirteen 
Ceo devices that were stable enough for detailed study 
(e.^.. Fig. 3c). This type of phenomenon, which has 
been measured previously in metal and semiconductor 
quantum dots j34l ES] , can be understood [20 in terms 
of excited-state-to-excited-state (ETE) transitions within 
the spectrum of energy levels in N^Ceo- The idea is il- 
lustrated by Fig. 4c, d. Here ground-state tunneling at 
low magnetic field corresponds to the transition (labeled 
A) from the q = -2, S = 3/2 state to g = -1, 5 = 1. For 
a sufficiently large applied voltage V, the tunneling of an 
electron back into the molecule can lead to a transition 
(labeled b) from the q = -1^ S = 1 state to the excited 
q = -2^ S = 3/2 multiplet, rather than the ground state 
6* = 3/2 multiplet. If this excited state does not relax 
before the next electron tunnels out of the molecule, this 
will produce an ETE conductance peak corresponding to 
a transition (B) between the excited ^^ = -2, 5' = 3/2 
multiplet and the excited q = -1^ S = 2 multiplet. The 
reason that transition B cannot extrapolate indefinitely 
to lower voltage is that the applied bias must provide 
sufficient energy to permit the transition b for transi- 
tion B to be possible. As support for this interpretation, 
note that the level diagram in Fig. 4c, d requires that 
Ec = Eb -\- {El) — Ea). Using the measured values of 
the excitation energies (see the caption for Fig. 3), this 
equality is satisfied within experimental uncertainty. 

Another successful feature of the calculation is that 
it can account for the magnetic-field dependence of the 
tunneling spectrum shown in Fig. 2a. In particular, note 
that this spectrum has two conductance peaks that as 
a function of magnetic field intersect the lowest-voltage 
transition but then terminate so that they do not pro- 
duce a level crossing - this happens near 1.5 T for tran- 
sition B and 4 T for transition C. Only near 7 T is 
there a true level crossing for transition D that corre- 
sponds to the spin-state transition discussed above. The 
"non-crossing" levels can be understood as another con- 
sequence of nonequihbrium ETE transitions. For volt- 
ages less than the ground-state to ground-state transi- 
tion, the first electron tunneling event that initiates cur- 
rent fiow is not allowed, so that the excited states cannot 
be populated. 

One difference between the calculated spectrum in Fig. 
4a, b and the measurements in Fig. 2a and 3a is that the 
calculated spectrum contains two more lines at positive 
bias (7 lines at 5 = instead of 5). The extra lines (D' 
and E') correspond to non-equilibrium ETE transitions 



from the excited q = -2, S = 3/2 state to the two g = -1, 
5* = 2 states (i.e., they are lower-energy satellites of tran- 
sitions D and E) . We suspect that these extra transitions 
are present but unresolved in the experimental spectra 
because the linewidths of the experimental transitions 
increase as a function of energy and the ETE transitions 
have weaker amplitudes than transitions originating from 
the ground state, thereby making them difficult to distin- 
guish at higher energies. The level broadening with in- 
creasing energy is not included in our simple model, but 
if we add in the measured broadening by hand, the tran- 
sitions is completely obscured and E' becomes only 
marginally resolvable. 

The capability to make electrical contact to single mag- 
netic molecules without destroying their magnetic char- 
acter provides the foundation for more advanced strate- 
gies to manipulate spin states. We have shown that the 
charge and spin values of the molecule can be identi- 
fied, that spin excitation energies can be measured, and 
that the tunneling spectrum of N^Ceo can be under- 
stood qualitatively within an elementary model. The 
next challenges will include developing the means to co- 
herently populate desired superpositions of spin states 
(using, for instance, applied microwaves O |36]), devel- 
oping accurate techniques to read out spin information, 
and achieving improved control over the interaction be- 
tween the magnetic molecule and its environment so as 
to minimize device-to-device variations. 



METHODS 

Synthesis: N^Ceo was synthesized by continuous 
nitrogen ion implantation into freshly sublimed fullerene 
layers [26 with a yield (N@C6o:C6o ratio) of -IQ-"^. 
The N@C6o contained in the harvested product was 
enriched and purified by multi-step high-pressure liquid 
chromatography (HPLC) [37]. Purity was checked by 
UV-Vis absorption, HPLC, and electron spin resonance 
and found to be better than 99.5%. 

Electrochemistry: We have confirmed, using electro- 
chemical measurements, that the interaction in N^Ceo 
between the N atom and the extra electrons on Ceo is 
relatively weak. Cyclic voltammograms were performed 
on N@C6o or Ceo films formed by drop casting from 
toluene solutions onto a 3 mm Pt disc electrode, polished 
to a mirror finish, used as the working electrode. A 
Ag/AgCl (saturated NaCl) and a Pt coil were used as 
the reference and counter electrodes, respectively. The 
scans, carried out at a sweep rate of 500 mV/s, were done 
in 0.1 M tetrabutylammonium hexafiuorophosphate in 
acetonitrile in a single-compartment cell. The solution 
was deoxygenated via N2 bubbling for 10 min before the 
sweeps, which were run in a N2 atmosphere. We found 
that the first two reductions of N^Ceo occur within 25 
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mV of those of Ceo- The 25 mV shifts are not related to 
the nature of the molecules and can be attributed to the 
varying rates of dissolution of the reduced species due to 
slight variations in the film morphology. This arises as a 
result of the fact that while the neutral forms of Ceo and 
N@Ceo are insoluble in acetonitrile, the reduced forms 
are not. 

Sample Fabrication: Our single-molecule transistors 
were fabricated following aprocedure similar to those 
used previously [H [5] . First we made an Al gate electrode 
16 nm thick and 2 fim wide on an oxidized Si wafer, and 
exposed the Al to air to form a thin insulating oxide. 
On top of the gate electrode, we fabricated continuous 
Ft wires with widths of approximately 150 nm and 
thicknesses of 10 nm. The chips were cleaned with an 
oxygen plasma and immediately covered with 25 /iL of 
either a 0.1 mM solution of N@Ceo for 2.5 minutes or 
(for one type of control sample) a 0.5 mM solution of 
Ceo ill toluene for 1 minute. Then the excess solution 
was blown off the chip, and the deposition process was 
repeated. We found that repeating the deposition pro- 
cess produced a more convenient yield of single-molecule 
devices. After the molecules were deposited, we cooled 
to cryogenic temperatures and broke the wires using 
electromigration in a circuit with small series resistance 
[17], thereby forming nm-scale gaps in which a molecule 
was sometimes trapped. The breaking voltages for our 
Ft wires were in the range 1-1.5 V. 

Low-temperature measurements: The mea- 
surements were performed in a dilution refrigerator 
equipped with room-temperature low-pass filters and 
low-temperature copper powder filters on all electrical 
lines. The electron temperature calibrated on quantum 
dot samples is < 100 mK. The DC current was measured 
as a function of swept source- drain voltage at fixed gate 
voltage and magnetic field (without a lock-in amplifier), 
and then the value of the gate voltage or magnetic 
field was stepped. The differential conductance was 
calculated by numerical differentiation. 

Theory: The differential conductance dl/dV is calcu- 
lated using the rate-equation approach in the sequential- 
tunneling approximation [20 . The electronic molecular 
Hamiltonian reads 

Hei = {e-eV*)^ala^^Ua\a^a\ai-Jse-SN-B{sl^S'N) 

where e is the on-site energy of electrons with spin cr =t, i 
created by aj., eVg is the electric potential at the posi- 
tion of the molecule, U is the Coulomb repulsion, J is the 
exchange interaction between the nitrogen spin Sat and 
the spin Sg = Xlacr' ^ti^cr^cr' /'^)^cr' clcctrous in the 

LUMO, and B is the magnetic field with a factor g/j^B 
absorbed. Only one LUMO is included, the other two 



are assumed to be split off due to molecular deformation 
and to not to affect dl/dV. The asymmetric contact ca- 
pacitances and the fact that the bias voltage is applied 
to the (arbitrarily labeled) left electrode are incorporated 
by writing Vg = aVg-\-/3LV with the applied gate voltage 
Vg and bias voltage V, where a = 0.15 and /3l = 0.25 
give the best fit. Asymmetric tunneling amplitudes Il^r 
between the LUMO and the left and right electrodes are 
assumed; a ratio tL/tR = 0.4 is used. In addition, we 
include low-energy excitations of energy hoOn in the n- 
charge state, which in this description for the sake of 
definiteness we assume to be vibrational. The transition 
rates between vibrational states contain Franck-Condon 
matrix elements [38 Mnrn,n'm' (with electron numbers 
n, n' and oscillator quantum numbers m, m') with rel- 
ative values of the nonvanishing components M2o,io = 
M2o,ii M2i,ii M22,io 1 and M2i,io 0.7; the 
matrix is symmetric. The other parameters used in the 
calculations are e = -3.23415, U = 3.1, J = -0.0004, and 
ksT = 0.00002 ah in eV. 
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